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ABSTRACT 


An  analytical  and  experimental  program  was  conducted  to  develop  acoustic  fatigue 
design  crireria  for  aircraft  structures  subjected  to  intense  noise  in  a  high  temperature 
environment.  Equations  for  the  dynamic  response  of  a  buckled  panel  were  formulated 
for  simply  supported  boundary  conditions  using  large  deflection  plate  theory.  Random 
amplitude  acoustic  fatigue  testing  of  representative  aircraft  structure  was  accomplished 
at  temperatures  up  to  60(PF  to  provide  data  for  correlation  with  the  analytical  results. 
Empirical  design  criteria  were  developed  in  the  form  of  equations  and  nomographs  for 
predicting  the  thermal  and  dynamic  response  of  aircraft  structures  subjected  to  combined 
environments.  The  empirical  design  criteria  are  presented  in  handbook  format  for  design 
use;  examples  and  computer  programs  are  also  presented. 
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I  -  INTRODUCTION 


A  number  of  investigations  have  been  conducted  in  recent  years  to  develop  and  refine 
acoustic  torigue  design  criteriu  for  aircraft  structures,  S c  o  cf  these  I" vestigotions  hove 
refined  the  design  criteria  for  s . if f ened-sk in  struf  .7.4'-s,  or.d  have  helped  to  reduce  the 
uncertainty  involved  in  predicting  aynamir  response  and  fatigue  c hoiac ter ist ic s  of  the:  ; 
structures  wnen  subjected  to  propulsion  system  or  oerodynumiv  noise  ut  ambient  temperatures. 
However,  when  unusual  structural  configurations  or  en vbonmon to i  'onc'itions  are 
encountered,  the  existing  design  methods  become  less  applicable  and  the  judgment  of  .‘he 
design  engineer  must  be  relied  upon  more  heavily.  Structural  applications  are  commonly 
encountered  reday  i  i  the  neai  noise  fic-ld  of  an  operating  turbojet  cr  turbeten  cnaine  where 
severe  noise,  liigti  temperatures,  sta'ic  loading,  and  vibrator'  uttet  occur  simultaneously 
or  in  conjunction  with  each  other. 

This  program  extends  the  basic  design  technology  for  ambient  eniper  gtu:  c  aluminum 
structures  to  include  the  effects  of  simultaneous  application  o.  thermal  and  acoustic 
environments .  A  comp! ete  and  derailed  desoi i ption  of  ibis  progi am  i s  presented  . n  AFFDL- 
T R— z7 3—  1  55,  PART  I,  including  on  anclyrical  development  for  the  dynamic  response  of  he  *ed 
structures  before  and  after  theimal  buckling.  The  primary  purpose  of  the  analytical  effort  v.n;, 
to  identify  tire  parameters  which  describe  tire  structural  'e.nonse;  fher  'he  data  re'  la¬ 
ments  of  the  exCu-imentol  p.ogiam  we  e  dei  fied  ui  detail .  Measured  data  were  correlated 
with  ti  e  analytical  .esults  to  estcblisli  empirical  desigr*  critf-'a  in  the  form  of  equations, 
nomc„  aphs,  and  computer  programs. 

Th^-  l  jnificant  results  cf  this  investigation  are  abstracted  in  this  document  as  an  aid  to  ‘he 
arr^  crtt  des’gne'  .  The  design  met.  odoiogy  is  presented  in  a  convenient,  abbreviated  format 
to  simplify  routine  use.  Tlv  v  governing  assumptions  ond  conditions  are  discussed  in  Section 
li. 

lne  oho i ng  sections  dtirneofe  three  av‘  table  methods  of  application: 

Section  III  -  Design  Relations  for  Manual  Calculation 
Secrion  IV  -  Des'gn  Nomographs  for  Graphical  Calculation 
Section  V  -  Computer  Programs 

The  limiting  conditions  to  the  design  methods  are  discussed  in  section  VI  to  give  the  de'.gn 
engineer  a  quick  guide  as  to  the  applicability  of  these  methods  to  a  particular  design  . 
Appendix  I  summarizes  substructure  sect.un  properties  as  on  oid  in  the  substructure  design  . 
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II  -  ASSUMPTIONS  AND  CONDITIONS 


The  ana lyrical  development  of  AFF DL -T K -/3- 1  5S,  Port  I,  is  based  on  certain  simplifying 
assumptions.  Many  of  these  initial  assumptions  were  negated  by  the  development  of  the 
empirical  relations.  The  significant  remaining  assumptions  ore  that: 

o  The  temperature  over  the  surface  ol  o  panel  bay  is  uniform. 

a  The  substructure  is  temperature-independent. 

o  The  structural  exciting  force  is  i  andom  amplitude  broad -band  noise  with  a 
Gaussian  distribution  ot  ampli  tudes. 

o  Dominant  structural  response  occur1  only  in  the  fundamental  mode  of  the  panel. 

The  first  of  these  must  be  closely  approximated  in  practice,  or  unsymmetric  buckling  modes 
will  occur  and  the  empirical  relations  will  not  be  val id .  The  second  condition  implies 
shoit  time  duration  and  localized  exposure  tc  the  thermal  environment,  so  that  the  sub¬ 
structure  remains  at  a  considerably  lower  temperature  than  the  skin.  The  third  condition 
implies  thet  the  design  criteria  are  applicable  to  oircruft  po//ered  by  jet  engines,  and  ore 
net  app  icohle  to  aircraft  powered  by  reciprocating  or  turboprop  engines.  The  last  of 
these  conditions  restricts  the  designs  to  consideration  of  single  mode  response  which  can 
prove  erroneous  for  structures  which  exhibit  significant  multi-mode  response. 

The  panel  size,  nomenclature,  and  sign  convention  for  the  design  equations  are  shown  in 
Figure  1  .  This  i  del  represents  a  single  bay  of  a  multi-bay  stiffenea-skin  structure. 

Figure  2  is  representative  of  a  continuous  multi-bay  structure,  where  the  model  iri  Figure  I 
represents  the  center  bay.  Throughout  this  handbook,  the  simplification  a  =-  02,  b  -  b2 
will  be  used  to  oescribe  the  center  bay  size. 

The  center  bay  (  r  bay  of  interest)  is  shown  in  Figure  3  i'olated  from  the  remainder  of  the 
structure.  The  height  of  the  stiffening  member  olong  the  long  side  of  the  panel  bay  is 
defined  as  d.  Throughout  this  reoort  the  panel  bay  aspect  ratio  is  assumed  to  be  greater 
than  one  (b/a  £  1 ) . 


UNIFORM  TEMPERATURE,  T 


b)  SUPPORT  STRUCTURE  GEOMETRY 


FIGURE  2.  MULTI-BAY  FLAT  STIFFENED  PANEL  CONFIGURATION 


FIGURE  3,  STIFFENED  PANEL  BAY 


II!  -  DESIGN  RELATIONS 


Empiricol  design  relofions  for  ombient  and  elevated  temperatures  are  presented  in  this 
section.  Ambient  temperature  design  criteria  for  aluminum  structures  are  from  AFFDL- 
TR-67-1562  and  AFFDl-TR-71 -I073;  elevated  temperature  design  criteria  are  from 
AFFDL-TP.-73-I55,  PART  I  .  '  The  source  of  design  equations  and  information,  other  than 

Reference  1,  will  be  acknowledged  in  this  and  the  following  sections  by  a  reference  in 
parenthesis.  No  ottempt  is  made  to  cross-reference  equation  or  figure  numbers.  All 
unreferenced  criteria  are  from  AFFDL-TR-73-! 55,  PART  I  . 


A.  Anb’ent  Temperature  Design  Criteria 


The  design  criteria  for  ambient  temperatures  are  unchanged  from  existing  criteria.  Only 
the  dynamic  response  of  the  structure  is  involved  in  the  design  so  long  as  the  ambient 
tcmperoture  state  does  not  cause  buckling  of  the  skin. 


1  .  Skin  Design 

The  skin  design  criteria  of  AFFDL-TR-67-1 56  are  volid  for  aluminum  structures  at  ambient 
temperatures.  The  dynamic  stresses  at  the  midpoint  of  the  long  side  will  be  greater  than 
those  at  the  midpoint  of  the  short  side  if  the  panel  responds  in  the  fundamental  mode.  The 
dynamic  stress  at  this  point  is  (Reference  2): 


o  -  1,62  x  10 

A 


-sTeI'^o1'25  5(f)  (b/o)1'75 

IvJ  ^T7J  •t>', 


Dynamic  stresses  at  this  point  are  also  given  by 
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The  fundamental  mode  frequency  for  a  single  bay  of  a  multi-bay  structure  is  given  by 
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The  aspect  ratio  parameters  are  defined  as 


F  ,  j  -  b/a  +  a/b 


(4a) 


AR  =  3(b/a)2  +  3(a/b)2  *  2  (4b) 

while  the  ocoustic  pressure  density  of  the  fundamental  mode  is 

L(S  /20)-9] 

♦(0  =  2.91x10  -  psi/v  Hz  (5) 

Comparison  of  these  equations  shows  that,  for  identical  configurations,  Equation  (1)  gives 
higher  stresses  than  (2)  .  Since  both  ore  based  on  the  same  analytical  model,  the  difference 
is  in  the  data  on  which  these  empirical  relations  are  based.  These  empirical  relations  can 
be  considered  as  bounds  for  predicting  dynamic  stresses.  The  form  of  the  latter  equation 
lends  itself  to  much  easier  solution. 


Dynamic  stresses  calculated  by  these  equations  are  used  with  suitable  random  amplitude 
fatigue  curves  to  estimate  fatigue  life.  Figure  4  contains  rondom  amplitude,  reversed 
bending  fatigue  curves  for  7075-T6  aluminum  and  6AI-4V  annealed  titanium  at  ambient 
temperature.  These  Curves  v/ere  developed  from  coupon  fatigue  tests,  ond  the  riveted 
fatigue  data  are  representative  of  typical  riveted,  stiffened-s'kin  structure. 

EXAMPLE:  A  flat  aluminum  alloy  stiffened  structure  is  to  be  designed  to  withstand  an 
estimated^ sound  pressure  spectrum  level  of  120  dB  for  5  x  lO^cyclesat  ambient  temperature. 
The  skin  design  is  determined  by  the  procedure  described  below. 

Assume: 

o  Damping  intio:  £  =  0.012 
o  Panel  width:  o  =  4.75  inches 
o  Aspect  ratio:  b/a  =1.5  (AR  =  10.08) 


From  Figure  4,  the  dynamic  stress  for  riveted  aluminum  (mean  value  fatigue  curve)  cor¬ 
responding  to  a  design  life  of  5  x  10^  cycles  is  3.8  ksi  rms.  Using  Equation  (1),  the 
minimum  skin  gage  is  found  by 
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From  Equation  (3),  the  fundamental  mode  frequency  is 
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o)  ALUMINUM  ALLOY  7075-T6 
REVERSED  BENDING  ZERO  MEAN  STRESS 


b)  TITANIUM  ALLOY  6AI-4V  ANNEALED 


FIGURE  4.  RANDOM  AMPLITUDE  LOADING  FATIGUE  CURVES 
AMBIENT  TEMPERATURE 
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Substituting  this  into  Equotion  (2)  gives 
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0.013  in. 


The  thicker  of  the  two  skin  thicknesses  should  be  selected  if  skin  failure  is  to  be 
completely  avoided . 
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2.  Stiffener  Flange  Design 

The  acoustic  loading  on  me  surface  of  a  stiffened  panel  is  transferred  to  the  substructure 
predominately  by  a  transverse  shear  loading,  causing  the  open-section  stiffeners  to  bend 
and  twist.  The  stiffener  loading  is  reacted  along  the  skin -stiffener  attachment  (rivet)  line 
and  at  the  clip  attachments  to  the  frames.  The  resulting  flange  stress  is  given  by 
(Reference  3), 


v°-9 

where 

XX  22  X 2  22 

This  relation  is  valid  only  for  fhe  fundamental  mode  of  the  panel. 

The  above  flange  stress  is  used  in  conjunction  with  o  fatigue  curve  developed  for  flange 
failures  and  presented  in  Figure  5  (Reference  3). 

EXAMPLE:  A  flat  aluminum  alloy  stiffened  structure  is  to  be  designed  to  withstand  on 
estimated  sound  pressure  spectrum  level  of  120  dB  for  5x10^  cycles  at  ambient  tempera¬ 
ture.  The  stiffener  design  is  determined  by  the  procedure  described  below. 
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Assume; 


o  Damping  ratio:  C  =  0.012 
o  Panel  width:  a  =  4.75  inches 
o  Aspect  ratio:  b/a  =  1  .5  (F  ^  =2.17) 

from  the  previous  example,  the  skin  thickness  is  token  os  h  =  0.025  inch,  for  which  the 
fundamental  mode  frequency  is 

f  =  1  .052  x  )04  x  0.025 
o 

=  2od  Hz 


The  stiffener  is  a  zee  section  0.040  inch  thick  with  a  flange  width  of  0.75  inch  ond  height 
ot  i  .25  inch.  The  section  properties  give  a  value  of  1*  =  0.01255  inch^.  (See  Appendix 
I  for  section  property  calculations .) 


From  Equation  (6),  the  flange  stress  is 
~  =  0.9 


0.0121(7. 125)3(1 ,25)(2.9  x  10-3) 

'263 

1/2 

(.01255X2.17) 

.012 

1/5 


=  2.19  ksi 


rms 


From  Figure  5,  the  life  is  estimated  to  be  N  =  1  .0  x  10  cycies,  or  slightly  less  thjn  the 
design  requirement. 


The  above  procedure  should  then  be  repeated  using  o  thicker  or  deeper  zee  stiffener  until 
the  desired  life  is  achieved.  It  is  also  possible  to  reduce  the  stiffener  spacing,  thereby 
increasing  the  fundamental  frequency  and  perhaps  altering  the  excitation  (dependent  on 
the  excitation  spectrum  shape).  However,  it  must  be  reme.ibered  tbet  changing  the 
spacing  or  skin  thickness  alters  the  skin  stress,  resulting  in  a  conservative  skin  design. 


3.  Damping  of  Stiffened-5kin  Structures 

The  estimation  of  a  damping  value  for  structural  design  is  normally  based  on  previous  test 
data  for  similar  structure  and  the  experience  and  judgment  of  the  design  engineer.  To  aid 
the  individual  designer  in  this  task,  the  damping  data  from  this  and  other  programs^'^  were 
plotted  versus  frequency,  and  the  resulting  plot  is  presented  in  Figure  6.  All  of  the  da*a, 
except  those  from  Reference  4,  are  for  fundamental  mode  response  of  stiffened-skin 
structure.  Although  the  data  scatter  is  such  that  no  definitive  estimation  criteria  are 
derived,  some  general  observations  may  be  drawn: 


*  '  : 

-3 


4  A 
% 


s 

* 
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. . *m”*p*??*m . . . , . . . . . . 


o  ALUMINUM  PANEL  DATA 
□  TITANIUM  PANEL  DATA 


A  ALUMINUM  PANEL  DATA  -  REF.  2 
<>  ALUMINUM  BOX  STRUCTURE  DATA  -  REF.  3 
A  ALUMINUM  PANEL  DATA  -  REF  .  3 


o  The  data  of  Reference  4,  which  included  response  of  a  single  panel  boy  through 
the  (3,3)  mode,  were  modeled  by  an  equivalent  viscous  damping,  whereby  the 
damping  ratio  is  inversely  proportional  to  the  forcing  frequency.  The  data  were 
thus  represented  by  an  equation  of  the  form 


c  =  11.0/f 


1.13 


A  similar  curve  was  established  by  Reference  5; 

C  =  3.6/f 

Both  of  these  representations  ore  shown  in  the  curve  ol  Figure  6. 


o  A  similar  approximation  was  established  through  the  data  points  of  Figure  6,  with 
the  slope  taken  os  an  average  of  the  Reference  4  and  5  slopes.  The  equation  of 
this  approximation  is 


c 


3.5/f 


1  .07 


o 


The  mean  value,  for  the  fundamental  mode,  of  all  damping  data  shown  in  Figure  6 
is 


c 


mean 


0.016 


This  mean  value  can  be  used  in  lieu  of  definitive  test  data  for  a  particular  con¬ 
figuration  . 


o  The  upper  and  lower  bounds  of  all  measured  data  are  Cu  =  0.038  and  C  |  -  0.0056. 
Application  of  these  bounds  on  damping  ratio  will  provide  limiting  values  for 
dynamic  stresses. 

The  foregoing  may  be  useful  for  rough  estimates  of  structural  damping.  The  actual  choice 
of  damping  ratio  is  left  to  the  judgment  of  the  design  engineer,  with  the  above  offered  as 
a  guide . 


B.  Elevated  Temperature  Design  Ci  iteric 


The  elevoted  temperature  design  criteria  are  used  in  essentially  the  same  sequence  in 
which  they  are  discussed  in  the  following  subsections.  Since  more  thon  one  method  of 
application  is  available,  all  the  criteria  will  be  summarized  and  then  followed  by 
exam, pies  of  their  use  . 
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1  .  Skin  Buckling  Temperature 

The  critical  buckling  temperature  of  o  single  panel,  such  as  that  shown  in  Figure  3,  is 
defined  as  the  temperature  increase  at  which  the  skin  buckles  and  is  given  by 


5-25h  F!l 

I  =  — t ■; — t —  -  F  above  ambient 

c  a  ab(l  +  v) 


(8) 


The  temperature  ratio  is  then  defined  as 


r  -  T/l 


(9) 


where  T  is  the  temperature  rise  of  the  structure  above  ambient. 


M 


2.  Skin  Buckling  Amplitude 

The  maximum  skin  buckling  amplitude,  at  the  midpoint  of  the  panel  bay  shown  in  Figure  3, 
is  given  by 


f 

l  ‘ 


t 

I- 


W  -  2 . 50  h  F ! ;  75  1  r-^- 
o  1 1  l  p 


1/2 


-  inches 


where  R  is  on  aspect  ratio  parameter  defined  as 


2^2 


R  =  3[(5  -  v  )F|  j  -  2(5  -t  v)()  -•.)] 


(10) 


(11) 


3.  Thermal  Stress 


Thermal  stresses  due  to  in-plone  expansion  and  skin  buckling  ore  given  by  the  following 
equations  for  the  midpoints  of  each  side: 


o  Midpoint  of  the  panel  bay  long  side  (y  =  b/2) 

,2 

U.OZ  t  VI 

1  -u" 


a  = 
x 


c  _  0.82  EW' 

E  a  T  ^  o 

T 


ab(l  -v)  La 
o  Midpoint  of  the  panel  bay  short  side  (x  =  a/2) 

,2 


-(2-vV  v°l  jlO-3 

bJ| 


a  - 

y 


c  T  1.66EW' 

EaT  t _ o 

1  -  v  ,  ,  2. 

ab(l  -  v  ) 


a  n  2.  b 

r-  (2  -  v  )  +  v  - 
b  a 
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-3 


-  kSl 


ksi 


(12) 


(13) 


fl 


Thermal  stresses  must  be  computed  at  the  midpoint  of  both  sides,  since  the  short  side  stress 
is  usua  I  ly  gi  eater  than  that  at  the  c  enter  of  the  long  side  .  This  is  opposi  te  lo  the  mogni  - 
fades  of  the  dynamic  stresses  at  the  two  locations.  The  above  relations  can  be  simplified  to 


3  +  a 

1  x. 


7  =  a  t  a 

y  <  ) 


where  the  respective  parameters  are  the  in-plane  expansion  and  buckling  stresses  as 
inferred  from  Equations  (1  2)  and  (13) . 


4.  Elevated  Temperature  Frequency  Response 

The  fundamental  mode  frequency  at  a  temperature  increase,  T,  isgiven  by  the  following 
equations  for  the  indicated  temperature  ranges: 


f(r)  -  f  [0.60  +  0.40(1  -  r)  '  J  -Hz  (0  «  r  «  1) 
o 


-  f  [0.60  +  0 . 44(r  -  1) 
o 


(r  *  1) 


5.  Dynamic  Stress 


Dynamic  stresses  at  any  temperature  can  be  computed  by  the  fol lowing  equations  for  the 
indicated  locations: 


o  Rivet  row  at  midpoint  of  long  side 


o  =  3.60  X  10 

x 


/b\2  *(f)  rf(r) 
\  h  /  AR  [  ' 


o  Rivet  row  at  midpoint  of  short  side 


a  =  13.0  x  10 
y 


fef 1/2 

\  h  /  AR  L  C  . 


The  elevated  temperature  panel  response  frequency,  f(r) ,  must  be  used  for  these  computa¬ 
tions.  The  stresses  at  both  locations  must  generally  be  calculored  for  elevated  temperature 
applications  because  of  the  interaction  of  the  thermal  and  dynamic  stresses. 

EXAMPLE •  A  flat  structure  is  to  be  designed  for  a  service  life  of  100  hours  at  a  sound 
pressure  spectrum  level  of  140  dB  and  a  service  temperature  of  500°F  .  Stainless  steel 
PH15-7Mo  is  selected  os  the  alloy  to  be  used  for  this  structure. 


Assume; 

o  Panel  width:  0-6  inches 
o  Aspect  ratio;  b/a  -  3.0  (F^  -  3.33) 
o  Skin  thickness:  h  -  0,050  inch 
o  Damping  ratio:  C-  0.016 
o  Ambient  temperature:  8Q°F 

A  fatigue  curve  for  the  selected  alloy  at  the  design  temperature  was  obtained  from  MIL- 
HDBK-5B^  tc  give  the  effects  of  mean  stress  on  fatigue  life.  (The  compressive  mean 
stresses  were  extrapolated.)  This  axial  loading,  constant  amplitude,  fatigue  curve  was 
converted  to  an  equivalent  random  amplitude  fatigue  curve  (Figure  7)  using  the  method  of 
Reference  7. 


From  Mll-HDBK- 5B,  the  alloy  properties  are 
V  --(.277/336)  7. 1  7  x  10~4  lb-sec2/in4 

c  -  6.1  x  10  ^  in/in/°F  {a-  500°F 
E  -  29.0  x  10°  osi  :  80°F 
E  --  26.97  x  10  psi  A  500° F 


The  skin  buckling  temperature  is  calculated  from  Equation  (8): 


T 

c 


5.25(  .050)2(3.33) 
(6.1  x  10"6)(6)(18)(1.32) 


-  50.3  F  above  ambient 


The  temperature  ratio  is  then 


r  -  420/50.3  8.35 

from  Equation  (9).  The  aspect  ratio  parameter  R  is  calculated  from  Equation  (II)  . as 
R  -  3[4. 8976(3. 33)2  -  2(5 .32)( . 68)3  =  141  .2 


The  buckling  amplitude  can  than  be  computed  from  Equation  (10); 


V' 

o 


2.50(.050)(3.33)1  ,75[7. 35/141 .2]  /?  -  0.235  inch 


This  is  the  maximum  displacement,  at  the  center  of  the  bay,  caused  by  the  50C°F  thermal 
environment.  The  nexi  step  is  the  computation  of  the  thermal  stresses,  usir  q  ecuations 
(14). 
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(26.97  x  I06)f6,l  x  10_6)(42C)  _ 


(0.82)(26.97  x  1 


(6)(l8)(.8976; 

A 


0.68 

o\.2'35)2 


=  -101  .5  ksi 


|^(3  •  0)<  1  .8976)  +  ^  j 


10  =  73.1  ksi 


(J.66)(26.97x  10  ;(.235)^  fl-8976  .  ,  „w.,  J  ,„-3  . 

Vb=  - WW^, - [-T0-^-32)(3-O)J  10  "40*6k11 


ond  then 

a  -101.5+  73.1  -  -28.4  ksi 
x 

a  »  -101  .5  +  40.6  =  -60.9  ksi 
Y 

The  ambient  temperature  fundamental  frequency  is  found  from  Equation  (3)  as 

x  Tl/2 


,  (0.79)(3.33)(  .050) 

f0  Tsyrrsr 


29.0  x  10 


=  259  Hz 


(7.17  x  10  )( .8976) 

Then  the  frequency  response  at  a  temperature  of  500°F  is  computed  from  Equation  (15;  as 
f(r)  «  f  >'0.60  +  0.44(7. 35)’ -  1  .79  f  Mi^Hz 

o  o 

The  acoustic  pressure  density  corresponding  to  a  spectrum  level  of  140  dB  is 
4(f)  =  2.91  x  | o'- (>40/20) -9]  t  2  9  x  ]Q~2 
using  Equation  (5) .  From  Equation  (4)  the  aspect  rotia  parameter  AR  is 

AR  -  3(3. 0)2  -  3/(3. 0)2  +  2  -  29.33 

This  then  allows  the  computation  of  the  dynamic  stresses  using  equations  (16)  ond  (17): 


n  -  3.6  x  10 
x 


-4 


18  ]2 

'2.9  x  10‘2" 

’464  ' 

.050 

L  .j 

29,33 

_.0l6_ 

1/2 


-  7.86  ksi 


rms 


13.0  x  10 


-4 


jl.TT 

. 050 J  L 


2.9  x  10 
”29.53 


464  T/2 

L”oT6  J  =  3-,5ks! 


rms 


Enter  the  fatigue  curve  of  Figure  7  (Kf  -  4)  with  the  dynamic  stress  ox  =  7.86  ksirrnj  ond 
thermo  I  mean  stress  ay  =  -28.4  ksi.  This  combi  na  t  ion  gives  a  life  of  opprox  imate  ly 
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8 

4.5  a  10  cycles.  Using  the  y-direction  stresses,  Oy  =3.15  ksirmi  and  Oy  -  -60.9  ksi, 
gives  o  Ijfe  greater  than  10^  cycles. 

At  c  frequency  of  464  Hz,  the  life  is  then  determined  by  using  the  minimum  cyclic  life, 

Q 

i .  f  _  N  _4.5xl0  _  L1 

L,fe "  jmw> '  ‘ 

which  is  greater  than  the  100  hour  design  requirement.  The  design  is  optimized  by  itera¬ 
tion  i  n  the  above  procedure  to  decrease  the  skin  gage  or  increase  the  stiffener  spacing 
such  that  the  predicted  life  is  equal  to  or  greater  than  100  hours. 
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IV  -  DESIGN  NOMOGRAPHS 


A  useful  tool  for  the  design  engineer  is  the  design  nomograph,  which  graphically  displays 
an  equation  for  rapid  solution.  The  empirical  relations  of  the  preceding  section  were 
formulated  into  such  nomographs  arvd  are  presented  in  this  section. 


A.  Ambient  Temperature  Design  Criterio 


1 ,  Skin  Design 

Design  criteria  for  stiffened -skin  aluminum  and  titanium  structures  are  displayed  in  Figure  8 
(References  I  and  2).  This  nomograph  is  applicable  to  7075-T6  aluminum  and  6AI-4V 
titanium. 

EXAMPLE:  A  flat  aluminum  alloy  stiffened  structure  is  to  be  designed  to  withstand  an 
estimated  sound  pressure  spectrum  level  of  120  dB  for  5  x  10®  cycles.  The  skin  design 
is  determined  by  the  procedure  described  below. 

Assume:  o  Damping  ratio:  £  =0.012 

o  Panel  width:  a  =  4.75  inches 

o  Aspect  ratio:  b/a=l.5 

Enter  the  nomograph.  Figure  8,  with  the  design  life,  and  follow  through  the  nomograph, 
as  indicated  by  the  arrows,  to  obtain  a  skin  thickness  h  =  0 .032  inch  . 

The  fundamental  frequency  is  estimated  from  Figure  9  (described  in  the  following  subsection) 
us  fQ  =  340  Hz.  At  this  frequency,  the  service  environment  spectrum  level  is  checked  with 
he  spectrum  level  used  above.  If  necessary,  an  iteration  is  made  to  obtain  agreement. 
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NOTE:  ALUMINUM  ALLOY 

DESIGN  NOMOGRAPH 
FROM  REFERENCE  2 


FIGURE  8.  STIFFENED  PANEL  SKIN  DESIGN  NOMOGRAPH 
AMBIENT  TEMPERATURE 


2.  Fundomeritol  Mode  Frequency 


lhe  fundamental  mode  frequency  for  a  single  bay  of  a  multi-bay  structure  is  given  by  the 
nomograph  of  Figure  9.  This  nomograph  was  developed  for  a  constant  value  of  Poisson's 
ratio,  v  =  0.32.  The  chart  was  simplified  by  taking  advantage  of  the  essentially  con¬ 
stant  ratio  of  E/y  for  most  aircraft  structural  alloys.  An  average  value  of  E/y  -  3.98  x  10^ 
in^/sec^  was  used;  this  is  an  average  of  the  ratios  for  aluminum,  titanium,  stainless  steel, 
and  Inconel  alloys . 

EXAMPLE:  The  frequency  of  the  structure  used  for  the  preceding  example  is  desired. 

Panel  width:  a  -  4.75  inches 
Aspect  ratio:  b/a  =  1  .5 
Skin  thickness:  h=  0.032  inch 

Follow  the  arrows  through  the  nomograph  to  obtain 

f  340  Hz 
a 


FIGURE  9.  AMBIENT  TEMPERATURE  FUNDAMENTAL  FREQUENCY  NOMOGRAPH 


B.  Elevoted  Temperoture  Design  Criteria 

The  examples  of  nomograph  use  shown  in  this  subsection  are  from  the  example  given  in 
Section  111 .  B . 

1  ■  Skin  Buckling  Temperature 

The  critical  buckling  temperature  of  a  single  bay  of  a  multi-bay  structure  is  given  by  the 
nomograph  of  Figure  10.  This  nomograph  is  based  on  a  constant  value  of  Poisson’s  ratio, 
v  =  0 . 32 . 

The  temperature  ratio  is  then  defined  os 

r  =  T/T 

c 

where  both  the  skin  temperature,  T,  and  the  critical  buckling  temperature,  Tc,  are 
measured  relative  to  the  ambient  temperature. 

EXAMPLE:  The  skin  critical  buckling  temperature  is  desired  for  the  structure  of  the 
example  in  Section  III  .B, 

Panel  width:  a  =  6  inches 
Aspect  Ratio:  b/a  -  3.0 
Skin  thickness:  h  -  0.050  inch 

Coefficient  of  thermal  expansion:  a  =  6. 1  x  lQ-^  irv/in/°F 
(PH15-7Mo  Stainless  Steel) 

Follow  the  arrows  Through  the  nomograph  to  obtain 

T  =  50^  above  ambient 
c 

then  the  temperoture  ratio  is 
r  =  420/50  =  8.4 


24 


BUCKLING  TEMPERATURE  NOMOGRAPH 


2.  Skin  Buckling  Amplitude 

The  sk  in  buckling  amplitude  is  given  by  the  nomograph  of  Figure  1 1 .  This  is  the  moximum 
buckle  amplitude  on  o  single  bay  of  a  multi-bay  structure;  it  occurs  at  the  midpoint  of  the 
bay.  The  nomograph  was  developed  for  a  constant  value  of  Poisson's  ratio,  v  =  0.32. 

The  buckling  amplitude  is  dependent  on  the  temperature  ratio  as  given  by  the  nomograph 
of  Figure  10. 

EXAMPLE:  The  skin  buckling  amplitude  is  desired  for  the  structure  of  the  example  in 
Section  III . B. 

Skin  temperature:  T  =  420°F  above  ambient 

Panel  width:  a  ~6  inches 

Aspect  ratio:  b/a  -  3.0 

Skin  thickness:  h  =0.050  inch 

From  the  preceding  example,  the  temperature  ratio  is  r  =  8.4.  Enter  the  nomograph  and 
follow  through  the  arrows  to  obtain  the  buckling  amplitude 

W  =0 .24  inch 
o 
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TEMPERATURE  -  F  ABOVE  AMB1EN  T 


3.  Thermal  Stress 

Thermal  stresses  due  to  in-plane  expansion  and  skin  buckling  are  given  by  Figures  12 
through  14.  The  thermal  stresse'  are  then  the  sum  of  these  component  stresses  in  the 
applicable  direction: 

o  Midpoint  of  pjnel  long  side,  in  x  -direction 

o  =  JT  +  o 
x  T  xb 

o  Midpoint  cf  panel  short  side,  in  y  -direction 

a  =  cT  +  a 

y  T  yb 

The  in-r.lane  expansion  stress  is  compressive  (-),  while  the  buckling  stresses  ov,  and  o,., 

*b  /b 

are  tensile  (+);  the  sign  must  be  maintained  throughout  the  computation. 


The  value  of  Poisson’s  ratio  has  been  held  constant  at  v -  0.32  for  development  of  these 
nomographs. 

The  $k>n  buckling  amplitude  from  Figure  1 1  is  required  for  solution  of  the  thermal  stresses. 

EXAMPLE:  The  thermal  stresses  are  desired  for  the  structure  of  the  example  in  Secrion  III  .6. 
Skin  temperature:  T  =  420°F  above  ambient 
Panel  width:  a  =  6  inches 
Aspect  ratio:  b/a  =  3.0 

Coefficient  of  thermal  expansion:  or  *6.1  x  10  ir\,/irv/0F  (S  500oF 
Modulus  of  Elasticity:  E  -  26.97  x  10^  psi  (S  500 °  f 

Enter  Figure  12  with  the  required  parameters,  and  follow  the  arrows  to  obtain 
c  ^  =  -100  ksi 

From  the  preceding  example,  the  skin  buckle  amplituae  Is  0.24  inch.  Enter  Figures  13 
and  14  with  the  buckle  amplitude  WQ  =  0.24,  and  follow  the  arrows  through  the  nomograph 
to  obtain 

a  -75  ksi  ;  a  =  40  ksi 


o  =  -100  +  75  =  -25  ksi 

x 

o  =  -100  +  40  -  -60  ksi 

X 
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FIGURE  13.  THERMAL  BUCKLING  STRESS  NOMOGRAPH 
X  -  DIRECTION  STRESS 


4.  Eli-voted  Temperature  Frequency  Response 


The  fundamental  mode  frequency  decreases  with  increasing  temperoture  un*‘i  the  skin 
buckles;  the  frequency  then  increases  with  further  increases  in  temperature.  The  nomo¬ 
graph  of  Figure  15  gives  the  elevated  temperature  fiequency  response  as  a  function  of 
the  ambient  temperature  frequency. 

EXAMPLE :  The  elevated  temperature  frequency  response  is  desired  for  the  structure  in 
the  example  of  Section  III. 6. 

Skin  temperature:  T  -  420°F  above  ambient 

Panel  width:  o  -  6  inches 

Aspect  ratio:  b/a  =  3.0 

Skin  thickness:  h  =0.050  inch 

From  the  preceding  examples  in  this  Section,  the  temperature  ratio  is  r  -  8.4.  From 
Figure  9,  the  fundamental  frequency  for  this  configuration  is  found  to  be 

f  =  260  Hz 

o 

From  the  nomograph  of  Figure  15,  the  frequency  ratio  is 

f(r)/f  =  1.8 
o 

from  which  the  frequency  response  at  the  temperature,  T,  is 

f(r)=  1.8  f  =  1.8  (260)  =  468  Hz 
o 


3.  Fatigue  Life 


* 


The  nomogroph  of  Figure  16  it  used  for  elevated  temperature  applications  for  7075-T6 
aluminum  a.  temperatures  of  150°  end  300°F  and  6AI-4V  titanium  at  temperatures  of 
400°  and  60r)nF,  7his  nomograph  includes  thermal  mean  stress  effects  in  the  data, 
thereby  negating  the  need  to  compute  these  parameters. 

EXAMPLE:  A  flat  aluminum  (7Q75-T6)  structure  is  to  be  designed  for  a  service  life  of 
100  hours  at  a  sound  pressure  spectrum  level  of  120  dB  and  a  service  temperature  of  300°F, 

Assume:  o  Aspect  ratio:  b/a  -  3.0  (F^  =  3.33) 

o  Damping  ratio:  Q  =  Q,G16 
o  A.mbient  temperature:  80°F 

As  a  first  step,  assume  a  frequency,  at  the  service  temperature,  of  300  Hz.  Then  the 
life  >n  cycles  is 

N  =  (300  Hz)(100  Hrs)(3600  sec/Hr  )  =  I  .03  x  108  cycles 


Enter  the  nomograph  cf  Figure  16  with  this  life  and  follow  through  the  parameters  to  the 
skin  thickness  chart .  Several  spacing/skin  thickness  ratios  are  now  possible,  all  of  which 
wi  1 1  meet  the  design  life  goal .  Assuming  a  skin  thickness  of  h  -  0.050  inch,  the  panel 
width  is  found  to  be  a  *  5.0  inches. 


At  this  point  a  structural  configuration  it  defined,  however,  the  assumed  frequency  must 
be  checked.  Compute  the  fundamental  mode  frequency  at  the  ombient  temperature  using 
Figure  9,  (fQ  =  390  Hz)  then  compute  the  skin  critical  buckling  temperature  using  Figure  10, 
where  or  -  13.5  x  10"®  in/in/6F  from  Figure  18;  (T£  =  33°F).  The  temperature  ratio 
is  then 

r  =  T/T  =220/33  =  6.7 
c 


Figure  15  then  yields  the  frequency  ratio  f(r)/fc  =1.65,  and  the  elevated  temperature 
frequency  is 

f(r)=1.65f  --1.65(390)  =  645  Hz 
o 

Since  this  frequency  is  greater  than  the  assumed  frequency,  the  anticipated  life  will  be 
less  than  the  design  goal,  hence  an  iteration  is  necessary.  The  above  procedure  is  repeated 
using  the  calculated  frequency  of  645  Hz  at  the  design  temperature.  One  or  more  itera¬ 
tions  may  be  necessary  to  obtain  agreement  between  the  initial  and  final  frequency  (or 
design  life) . 


Several  spacing/ skiri  thickness  ratio  combinations  may  be  derived  using  this  method  and 
the  weight  of  each  calculated  to  obtain  a  minimum  weight  design. 


This  method  may  also  be  used  for  structural  temperatures  different  than  those  indicated  on 
the  nomograph  by  assuming  a  linear  relationship  between  the  temperatures  shown  and 
i  nterpolati  ng . 


1 

1 
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6.  Applicotion  of  the  Design  Procedure 


At  least  two  alternative  methods  of  application  ore  possible  using  the  eleveted  temperature 
design  criteria.  The  first  of  these/  and  the  most  direct,  involves  the  use  of  the  life  design 
nomograph,  Figure  16.  The  alloys  and  structural  temperatures  must  match  those  for  which 
the  nomograph  was  developed.  An  example  of  the  use  of  this  nomograph  was  presented  in 
the  preceding  subsection. 


a.  Mean  Stress  Fatigue  Curve  -  An  alternate  design  method  involves  the  use  of  fatigue 
curves  where  the  mean  stress  effects  are  known  (see  Figure  7).  This  method  can  be  employed 
where  the  alloy  or  temperature  does  not  coincide  with  those  of  the  Figure  16  nomograph. 

The  following  example  is  identical  to  the  example  of  Section  III.B  and  to  the  examples  used 
on  each  of  the  design  nomographs  of  this  section.  It  is  presented  here  to  summarize  the 
application  of  the  design  chart  section. 


EXAMPLE:  A  flat  stiffened  structure  is  to  be  designed  for  a  service  life  of  100  hours  at  a 
sound  pressure  spectrum  level  of  140  dB  and  a  service  temperature  of  500 °F.  Stainless 
steel  PH15-7Mo  is  selected  as  the  alloy  to  be  used  for  this  structure. 


Assume: 


o  Panel  width:  a  =  6  inches 
o  Aspect  ratio:  b/a  -  2  .0  (Fj  |  =  3.33) 
o  Skin  thickness:  h  =  0.050  inch 
o  Damping  ratio:  Q-  0.016 
o  Ambient  temperature:  80°F 


Figure  7  shows  a  random  amplitude  fatigue  curve  with  varying  mean  stresses  for  the 
PH15-7Mo  alloy. 

From  MIL-HDBK-5B,  the  alloy  properties  are 

y  =  (.277/386)  =  7. 17  x  10"4  lb-sec2/in4 

a  =  6.1  x  10"6  in/in/°F  @  500°F 

E  =  29.0  x  106  psi  @  RT 
o  r 

E  =  26.97  x  106  psi  ©500°F 

The  skin  critical  buckling  temperature  is  found  from  Figure  10  as  Tc  =  50°F  above  ambient. 
Then  r  =420/50  =  8.4  and  the  buckling  amplitude  is  WQ  =0.24  inch  from  Figure  II  . 

The  thermal  stress  components  at  the  midpoints  of  the  two  sides  are  found  from  Figures  12 
through  14,  or 

Cj  -  -100  ksi 

a  ■  75  ksi  ;  o  “40  ksi 
xb  yb 
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Then  the  thermal  (compressive)  stresses  are 

a  =  aT  +  a  =  -100  +  75  =  —25  ksi 
x  I  xL 


o  =  oT  +  a  -  -100  +  40  =  -60  ksi 

y  T  yb 

The  ambient  temperature  fundamental  frequency  is  fQ  -  260  Hz  from  Figure  9.  Figure  15 
gives  a  frequency  ratio  of  1.8,  corresponding  to  a  temperature  ratio  of  8.4;  then 

f(r)  =  I  .8  f  =  1.8(260)  =  468  Hz 
o 

The  dynamic  stresses  are  then  computed  by  Equations  (16)  and  (17): 
o  Midpoint  of  panel  long  side 

o  =3.60xl0~4  18  l2  2.9  xl0“2'  468  ''/2  =  7.89  ksi 

x  [0.050  29.33  0.016  rms 

o  Midpoint  of  panel  short  side 


0  =13.0xl0"4r  6  1  2  2.9  x  10“2  T  468  11/<2  =3.17  ksi 

y  0.050  29.33  0.016  rms 

-2 

The  value  of  t(f)  =  2.9  x  10  is  the  acoustic  pressure  density  corresponding  to  140  dB 
(Equation  (5)),  while  AR  =  29.33  is  the  aspect  rotio  parameter  (Equotion  (4b)). 


Enter  the  fatigue  curve  of  Figure  7  (Kf  -  4)  with  the  dynamic  stress  cx  =  7.89  ksirms  and 
thermal  mean  stress  cx  =  -25  ksi  ,  This  combination  gives  a  iife  of  approximately  4.5  x  10 


cycles.  Using  the  y-direction  stresses,  3y  “3.17  ksi  rms  and  Oy  =  -60  ksi ,  gives  a  life 
greater  than  10^  cycles.  At  a  frequency  of  468  Hz,  the  life  is  then  determined  by  using 


the  minimum  cyclic  life, 

LIFE  =  N  =  4.5  x  10B  =  267  Hours 


which  is  greater  thGn  the  100  hour  design  requirement.  The  design  can  be  optimized  by 
iteration  on  the  above  procedure  to  decrease  the  skin  gage  or  increase  the  stiffener  spacing 
such  that  the  predicted  life  is  equal  to  or  greater  than  the  100  hour  design  life.  Alterna¬ 
tively,  the  above  design  may  be  used  and  the  additional  estimated  life  used  as  a  safety 
factor . 


A  comparison  of  the  results  achieved  here  (using  the  design  charts)  and  the  results  of 
Section  I II.  8  (using  the  design  equations)  shows  a  difference  of  only  2  hours  (less  than  1% 
error)  on  the  final  life.  Hence  the  accuracy  of  the  graphical  melhod  of  solution  is  com¬ 
patible  with  that  of  the  computational  method. 
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After  optimizing  the  above  design,  the  stiffener  flange  details  can  be  determined  using 
the  method  of  Section  III. A.  The  elevated  temperature  response  frequency  calculated 
above  (i.e.,  f(r)  -  4b8  Hz)  should  be  used  for  this  computation.  The  fatigue  data  used 
in  the  skin  design  (i  .e . ,  Figure  7  for  this  example)  may  be  used  for  the  stiffener  design 
in  lieu  of  a  fatigue  curve  developed  specifically  for  stiffeners  at  the  elevated  temperature. 


38 


V  -  COMPUTER  PROGRAMS 


A  digilal  computer  program  is  presented  here  for  calculating  the  dynamic  response  and 
life  of  a  multi-bay,  fiat  <tiffened-skin  structure  exposed  to  simultaneous  acoustic  and 
thermal  environments.  Five  sub-proarams  ore  required  for  the  dynamic,  analysis  program. 
These  programs  were  developed  for  the  Uni.ac  1 106  computer  using  Fottron  V;  however, 
the  programs  can  readily  be  adopted  to  any  digital  computer. 

A.  Analysis  Program 

The  input  data  formot  for  the  analysis  program  >s  shown  in  Table  I  ond  the  input 
parameters  are  defined  in  Table  II.  The  computer  program  is  listed  in  Table  III  while 
Table  IV  contains  a  sample  of  the  output  format. 

6.  Sub-programs  Required 

The  following  sub-programs  are  required  for  this  analysis  program: 

ETEMP  (T,  IFF) 

ALPHA  (T,  IFF) 

SN  (SDYN,  STEMP,  TEMP,  CTF,  IFF) 

CTEMP  (TCALP,  TC,  IFF) 

PROP  (OPT,  B,  H,  T,  A,  RJ,  GAMAT,  XIP) 

The  input  and  output  parameters  for  each  sub-program  are  given  in  the  following  sub¬ 
sections.  A  listing  of  each  of  the  sub-programs  is  contained  in  Tables  V  through  IX. 
Since  each  of  these  sub-programs  are  either  functions  or  subroutines,  the  input  and 
output  are  controlled  by  the  coiling  program. 

1 .  Sub-programs  ETEMP  and  ALPHA 

These  functions  compute  modulus  of  elasticity  and  coefficient  of  thermal  expansion, 
respectively,  for  aluminum  and  titanium  alloys  as  a  function  of  temperature.  The  basic 
alloy  properties  are  from  MIL-HDBK-5B,  and  are  shown  in  Figures  17  and  18.  The  input 
parameters  are: 


T  -  Input  Temperoture  at  which  elastic  modulus  or  coefficient 
of  thermal  expansion  desired  -  °F 

IFF  -  Alloy  code, 

=  1  Titanium  Alloy  (6A1-4V  annealed) 

=  2  Aluminum  Alloy  (7075-T6) 
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b)  TITANIUM  ALLOY  6AI-4V 


FIGURE  18.  TEMPERATURE  EFFECTS  ON  COEFFICIENT  OF  THERMAL  EXPANSION 
(FROM  REFERENCE  6) 


2.  Sub-program  SN 

This  subroutine  computes  the  fatigue  life  of  aluminum  and  titanium  alloys  as  a  function 
of  dynamic  ui.d  thermal  mean  stresses.  The  input  parameters  are: 

SDYN  -  Dynamic  Stress  -  ksi  rms 

STEMP  -  Thermal  (or  Mean)  Stress  -  ksi 

TEMP  -  Temperature  -  °F 

IFF  -  Alloy  Code, 

=  1  Titanium  alloy  (6A1-4V  annealed) 

=  2  Aluminum  Alloy  (7075-T6) 

The  output  to  the  calling  program  is: 

OTF  -  Life  in  cycles  to  failure 


3.  Sub-program  CTEMP 

This  subroutine  computes  the  skin  critical  buckling  temperature  for  aluminum  or  titanium 
structures.  The  inpur  parameters  ore: 

TCALP  -  Product  of  critical  buckling  temperature  and 

coefficient  of  thermal  expansion;  computed  using 
Equation  (8). 

IFF  -  Alloy  code, 

=  1  Titanium  Alloy  (6A1-4V  Annealed) 

=  2  Aluminum  Alloy  (7075-T6) 

The  single  output  parameter  to  the  calling  program  is: 

TC  -  Critical  buckling  temperature  -  °F  above  ambient. 


4.  Sub-program  PROP 

This  program  computes  stiffening  member  properties  such  os  area  and  moment  of  inertia. 
The  basic  relations  are  from  Reference  3  and  are  shown  in  Appendix  I.  Two  different 
sectional  shapes  are  available,  a  zee  or  a  channel  section,  with  the  parameters  described 
in  Appendix  I. 

The  input  parameters  are 

OPT  -  Option  code  to  select  sectional  shape, 

-  0  zee-section 

=  1  channel  section 

8  -  Flange  width  of  stiffening  member  -  inch 

H  -  Height  of  stiffening  member  -  inch 
T  -  Thickness  of  stiffening  member  -  inch 
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The  output  parameters  to  the  calling  program  arc: 

2 

A  -  Cross-sectional  area  -  in 

RJ  -  St.  Venant's  Torsion  Constant  -  in^ 

GAMAT  -  Warping  constant  for  thin-walled  open-section 

beam,  with  the  pole  taken  at  the  shear  center  -  in^1 

XIP  -  Polar  moment  of  inertia,  referenced  to  rotation 
about  the  attachment  point  -  in^ 
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,AvLE  I 


TABLE  -  II 


NCASE 

IFF 


OPTX 

BX 

HX 

TX 

OPTY 

BY 

HY 

TY 

A1 

A2 

81 

B2 

TS 

RHO 

RNU 

DAMP 

PSL 

T 


DYNAMIC  ANALYSIS  COMPUTER  PROGRAM 
INPUT  PARAMETER  DEFINITION 

Two-digit  identification  number 

Alloy  identification  code, 

=  1  Titanium  Alloy  (6A1-4V  Annealed) 

=2  Aluminum  Alloy  (7075-T6) 


Input  parameters  defi...ng  stiffening  member  parallel  to 
x -direction  -  see  Sub-program  PROP  for  definition 


Input  parameters  defining  stiffening  member  parallel  to 
y-direction  -  see  Sub-program  PROP  for  definition 


Panel  bay  dimensions  -  see  Figure  2 


Skin  thickness  -  inch 

Weight  density  of  skin  and  stiffening  member  alloy  -  Ib/ir^ 

Poisson's  ratio  for  structure  alloy 

Damping  ratio  for  structure 

Sound  pressure  spectrum  level  -  aB 

Structure  temperature  rise  -  °F  above  ambient 
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TABLE  ill 


COMPUTER  PROGRAM  FOR  ELEVATED  TEMPERATURE 
DYNAMIC  RESPONSE  OF  STIFFENED  STRUCTURE 


1 

2 

3 

4 
b 
o 
7 
H 

CJ 

10 
1 : 
12 

13 

14 
lb 

1<j 

17 
1  H 
10 

2'J 

21 

23 

24 
2b 

2n 

27 

27 

/If 

S' 
3-  i 


0  o 
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C  THIS  PROGRAM  CALCULATES  THE  DYNAMIC  RESPONSE  OF 

C  A  NINE-BAY  FLAT  STIFFENED  PANEL  EXPOSED  TO  A 

C  UNIFORM  ACOUSTIC  PRESSURE  AND  A  UNIFORM  TEMP- 

C  ERATURE  RISE.  ROOM  TEMPERATURE  IS  BO  DEGRESS  F. 

C 

C  T  IS  A  TEMPERATURE  RISE*  ABOVE  ROOM  TEMPERATURE 

C 

C  SUBPROGRAMS  REQUIRED!  ALPHA » T * IFi  ) »  ETEMP ( T  * IFt  )» 

C  SN ( SO YN • STEMP  #  T  #  C  TF  » IF!  )  .  CTEMP ( TCALP # TC *  IF-  )  i 

C  AND  PROP (OPT.B.H.T.A.ZJ.WC.PIP) 

C 

C  FUNCTION  DEFINATION 

C 

F (B*A)=B/A*A/B 

R (B. A. PR/ =3.* ( <  S.-PR* •  2 1  * ( B/A+A/B ) * ■ 2-2.* (b. +PR ) 

1  * ( 1 . -PR ) i 

20.  READ(5.301  INCASE.IFi 

PEAD(b.302)OPTX.8X.HX.TX 
READ (b»302 10PTY.BY.HY .TV 
READ  ( b »  30  1  )K>' 

READ ( b#  303 ) A  1  .A2.B2.U1 
READ  (b.  303  ITS.  KHO.RflU.  DAMP 
READ(b.303)PSL.T 

C  INPUT  DATA  FORMAT  STATEMENTS 

301  FORMAT (212) 

302  FORMAT ( I2.3PB. 4) 

3  U  3  FORMAT (4FB, 4) 

C  CALCULATE  SUBSTRUCTURE  PROPERTIES 

CALL  PROP ( OPTX »BX  »HX  »  TX . AX  *  v j  »  WCX  #  X  I  ) 

CAL'  PROP  (QPTY  *  13  Y  »HY  »TY  »AY  *  Y  J  *  t:C  Y  »Yl  ) 

H  =  TS 
GM=RHO 
PR  — *  'IIU 

c  calculate  st i f  cu : r  STIF-  IIES  AND  mas 

RX  1  =  0  .-'bub  * A  1  .  A  1  *  X J/ <  WCX ♦ ( 1 .  »MR  ) 

RX2  =  0  .  Ob.  In  ,*  A 2  *  A  2  •  X  J/  i  V.‘C  X  *  (  1  .  HJK  )  ■ 

RY1  -0.  Ob:  )h  *131  Mil  tYJ/(V;CY*  1  1  .  ►PR  i  ' 

RY2=0.tib0o  •♦U2*il2»YJ/(WCY*  !  1  .  t-PR  ) 

SKX=Y/CX*  !  1  .  URX2K’.*  (A2/A1  )  ►  (  1  .  ii;xl  )  ’  /A2 

SKY=wCY+  (  1  .  +  KY2*  .  *  (  82/111  )  ♦  (  1  .  iRY  1  i  /LI 2 

i  (3=1 1  *i  i  m  : 

ST  R  =  4  73. 74  00*  (  I  ,  -f|<  *Pl<  )  V  I  SK  X  WY-  Y  )  /  (M3*  A  2  'll.  ) 
STP=STi-'/40 
A 3= ( A  1 /A2 ) ♦  3 
H  3=  (151/112)  ♦  3 
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TABLE  111  (CONT'D) 


GM='<  /3Bo. 

SKUrii  . 2b*GM  »i  l  A  2  'U2 '  (  1  .  *2  .  +  A3  +  2.  +B3  +  4  .  *A3  i  H 3  ) 

1  *‘.*,Bt)'i>.*GH  '  <  X  I +A2  •  11  .  +  2.  +A3)  /  (02*1)2  ) 

2  +YlU12t  (  1  .  »2.*H3)/(A2+a;m 
CALCULATE  COVE';  SHE  T  ST  IF  NES  AND  MAS 

F2  -F  ( B  2  *  A  2  ) 

F21=F  <B2* A  1  ) 

F 1 2=F ( B 1  * A2  ) 

FI  =F  (  ci  1  •  A  1  ) 

F1S=F2  +  F2  ♦  .+  (A1/A2  )  +  F2J  +  F21 
1  +2.*(U1/M2)»F12*F12+4.*(A1/A2)*(bl/U2)*F1  +  F 1 

RO  M  TEMP  ST  IF:  1 1ES 

GKU  =  2  .(124  37  +ETEMP  (UO  .  .  IF  )  *H3  '  (F  1S+STH  )  /  (  <  1  ,-PH+PH  ) 
1  *A2  U12  ) 

RO  M  TEMP  FREQUENCY 
F  li  =  i  ■ .  1  64+ SORT  (  GK  (J/SKM  ) 

CALCULATE  RO  M  TEMP  MEAN  SOUAHE  STUES'  RESPONSE 
AR=3.+ (B2/A2! +  2  +  3.+ (A2/B2)*  2+2. 

CONVERT  OB  TO  PSI 
SPL=2.01 *1 U. ♦ ■ (PSL/2U .-0. ) 

CALCULATE  RO  M  TEMP  DYNAMIC  STRESS  AT  X=0»Y=B2/2 
SXU=n.  3A*+B2  +U2  +  S0.  T  (FU/DAMP  )  +  SPL/  (H+H  »  AR  ) 

CALCULATE  RO  m  TE'^u  DYNAMIC  STKE$‘  AT  X=A2/2*Y=U 
SY 0-  1  •  30  + A2  + A2  +SOKT  (F  U/DAMP  )  +  SPL/  (H*i  U  AR  ) 

CONVERT  STRES  FROM  PSI  TO  KSI 
SX!l  — SXu/l(! 

SY(l  =  SYli/l  i  '.ii 

CALCULATE  RO  M  TEMPERATURE  LIFE 
CAL  Si  I  (  SX  U  »  0  .  !J  #  BU  .  1 1  #  CTF  1  *  IF'  ) 

CALI  SN l S Y  0 ► R •  '  i  *  Bu • U  *  CTF2  *  IF1  ) 

X1=A2/-’. 

Y  1=0  .  ij 
X2=0  .  U 
y2=b  2/;: . 

STLMP  =  0 .  * i 

PRINT  HO  M  TEMPERATURE  RESPONSE 
W'R  I TE  (  b  *  4 Or  ) 

GO  T0(201 *202) * JF+ 

201  WRITE ( 6*  40b )  NCASE 
GO  TO  203 

202  WRITE lb*4uo>  NCASE 

203  WRI TE ( b  *  41 0 )  PS-  +  T  ' 

WRITE (b*4Ib) 

WR  I T  E  (  6  *  4  1 6  )  F  U 


TABLE  III  (CONT'D) 


1  ‘II 

ivH  I  TL  (  ci  »  42 0  ) 

■»1 

WHITE  to. 42b) 

u2 

WHITE (0.430 )  X2.Y2.SXU#  STEKP  » C  1  F  1 

“3 

WHITE  (6.4  30  )  XI  .  Y  1  #GYU  .NTLMP.CTF2 

Ml* 

C 

Ob 

C 

THERMAL  STHEb  EH  LCTS 

l)n 

C 

“7 

H  2  =H (02.A2.HH ) 

PH 

H21=H (H2.A1  .PH  ) 

L  1 

H 12=R ( B 1 . A2.PR  ) 

in 

HI  ' =R  <B1 .A1 »PR  ) 

Hi  ’ 

H2S=F2  ’*  (F2  +2.*  (  A1  XA2  )  *  •  2  i  F2 1  f  2  .  *  (  b  1  /i  J2  )  ♦  2  F  12 

i  ii;.j 

1  +4.*(A1/A,_j»  2MBl/o2)*  2* FI  ) 

1  ii  3 

HST  =R2  +2.*A3*R12+2.+B3'P21  f’4.+A3»B3  >R  1 

1  !'4 

HOC  (F  1S  +  '  TR  )  /F2S 

1  ub 

C 

CALCULATE  CRITICAL  TEMPERATURE  RISE.  TCR 

1  1)0 

T C  ALP-b .  2b*l  1  • :  l  -F2  /  (  A2HJ2  Ml.  fPR  )  , 

1  U  7 

CAL  CTEmR(TCALF>»TCR.IF.  ) 

1  OH 

C 

*  MOTE  *  TCA  Al  ID  HS  ARE  BUCKLING  TEMPERATURE 

100 

C 

AMO  TEMPERATURE  RATIO  FOR  AN  EOUAL 

1  ;  0 

C 

SIZE  SIMPLE  PANEL.  R9  IS  TEMP  RATIO 

1 

C 

FOR  NIUE-BAY  PANEL 

i  ■> 

A  C— 

TCA-TCR/RO 

1  3 

RS-T /TCA 

1  ■  4 

R9-T / T  CH 

1  '  b 

TACT=T+80.0 

1  '  u 

C 

CALCULATE  MATERIAL  PROPERTIES  AT  TEMPERATURE 

1  7 

ES=ETEMP (TACT . IFF  ) 

1  '  H 

ALP=ALPHA (TACT » IF:  ) 

1  '0 

D=0 . 083o*LS  *H3/  ( 1 ,-PR+PR  ) 

120 

C 

CALCULATE  response  frequency  at  temperature .t 

121 

SKT=D*F2S*R0/ <A2*B2 ! 

12 

FQT  =  0 . 8 09* SORT (SKT/SKM) 

123 

C 

♦iNOTE*.  FOT=FO.  ROOM  TEMP  FREQUENCY 

124 

STLlN  =  -E5*ALP*T/< 1 . -PR ) / 1 0 un  .  0 

12b 

IF (RS-R0120b.20b.210 

12o 

C 

PRE-BUCKLED  RESPONSE 

127 

20b 

FTEMP=FQT* ( 0 .fcO+O . 40+SQRT ( 1.-R9) ) 

128 

SXT=STLIN 

129 

SYT=STLIN 

130 

W 0=0.0 

131 

GO  TO  21b 

132 

C 

post-buckled  RESPONSE 

133 

210 

FTEMP=FQT* < 0 . 60+0 . 44*SQRT (R9-1 .  )  ) 

134 

C 

CALCULATE  PLATE  BUCKLING  AMPLITUDE.  WO 

13b 

W0  = ( 3.37-0 .20 *R0 ) *H*SQRT  <F2S*R0* (R9-1 . ) /RST  ) 

48 


TABLE  III  (CONT'D) 


CALCULATE  THi  iEHMAL  STRES' ,ES 
C1=1./(A2*B2*<1.-PR*PR) ) 

SXT-STLIN-*-0.81*ES*Cl*(  <2.-PR*PR  )  *62/A2+A2*PR/B2  ) 

1  *W0+W0/10Uw  a  0 

SYT=STLlN  +  1.36*F.S*Cl*(PR*B2/A2+t2.-PR*PR  >*A2/B2) 

1  *WO*W0/10UU . 0 

CALCULATE  DYNAMIC  STRESS 
CONTINUE 

C2=SQRT (FTEMP/FO) 

SXO=C2*SXO 

SY0=C2*SY0 

CALCULATE  ELEVATED  TEMPERATURE  LIFE 
CALL  SN ( SXO  »  SXT  »TACT  #  CTF 1 » I F h  ) 

CALI.  SN  <SY0#SYT  .TACT  #CTF2»  IF,  ) 

PRINT  ELEVATED  TEMPERATURE  RESPONSE 
WRITE (b»43b) 

TA=TCR  +  80.t) 

WR ITE ( b  »  4i  0  )  TCR 
WRITE(b#4.b)  WO 
WR ITE  (  b  »  4 16 )  FTEMP 
WRITE (b *420) 

WR  I TE ( b  *  42b ) 

WR I TE ( 6  *  430 )  X2#Y2#SXU#SXT»CTF] 

WRITE (b *430 )  XI i Y1 #SYO#SYT »CTF2 
GO  TO  2  U '  ■ 

FORMAT  STATEMENTS  FOR  OUTPUT  DATA 
FORMAT (* 1 ' #2bX. 'DYNAMIC  RESPONSE  OF  A'»/#19x» 

1  'NINE  —NAY  STIF'  ENI  D  PANEL  EXPOSED  T0'./*21X# 

2  '  ACOUST I C  EXCITATION  AMU  HEATING'#/) 

FORMAT  <29X  »  'DAT'.  C  ASE  *  »  I  4  #  /  •' »  2  7X  #  'MATERIAL  *.  TITANIUM') 
FORMAT  (  29X  » 'DAT  CmSE  '  #  I  4  #  /  /  #  2  /  X  #  'MATERIAL  :  ALUMIIIU.,,1 
FORMAT (bX» 'EXCITATION  SPECTRUM  LEVEL  -  ' ♦ F4 . U » 1 X . « DB » * 
13X. 'TEMPERATURE  INCREASE  =  '  -F4 . U » IX #  'DEG  .  F»./<) 

FORMAT ( 24X . 'RO  M  TEMPERATURE  RESPONSE'#/) 

FORMAT  (  20X  »'  FI  li  lOAUF'  :T  .1  FREONCNCY  ='»F7.1#»  UZ  »  .  /  ) 
FORMAT  (bX# 'STRES  AT  PO I  NT  '  #  3X  .  »  D  YIIAM I C  STRES  '»3X» 

1 'THERMAL  STRES  *  #  3X  # 'CY'.LES  TO  FAILURE') 

FORMAT  <HX#  '*•’  #7x  .  '  y  •  .  i  X.'KSI'.14X#*KSI*#/) 

FOR'- AT  (  NX  #  F  b  .  2  »  3X  #  E'  S  «  2  #  Sx  » Ffl .  3  #  UX  #F  H  .  3  •  1  li  X  •  1 1'E ,  2  #  /  ) 
FORMAT  (/  #2  y.'LL  VATED  TEMPERATURE  M'CSPO'  ISE  »  #  /  I 

FORMAT  I 10X # 'DUCKL I HO  TEMPT  NATURE  ='.FH,2#'  DEG.  F» 

].'  AHOVE  RO  Mi  TE-iPE.RAT'  ii  :E  •  »  /  ) 

FORMAT U Hx . 'BUCKLING  AMPLITUDE  ='»FH,4# 

1  '  INCHES' ./ ) 

El  iU 


4? 


TABLE  IV 


OUTPUT  FORMAT  FOR  DYNAMIC  ANALYSIS  COMPUTER  PROGRAM 


DYNAMIC  RESPONSE  OF  A 
NINE-BAY  SI  IP i  EN'  D  PANEL  EXPOSED  TO 
ACOUSTIC  EXCITATION  AND  HEATING 

DATA  CASE  4 

MATERIAL  :  ALUMINUM 

EXCITATION  SPECTRUM  LEVEL  =  135.  DB  TEMPERATURE  INCREASE  =  £0.  DEG.  F 


RO< M  TEMPERATURE  RESPONSE 
FUNDAMENTAL  FRCQU'.NCY  =  175.  -  HZ 


STRES 

X 

AT  POINT 

Y 

DYNAMIC  STRES 
KSI 

THERMAL  STRES 
KSI 

CY-  LES  TO  r  /\  1  LURE 

.00 

6.  u: 

6  .<J  76 

.  u 

6.43+1 u 

3.0 

.  U  i: 

5.485 

.U 

1 .u3+i7 

ELI VATED  TEMPERATURE  RESPONSE 

BUCKLING  TEMPERATURE  =  1 '■ ,0>;  DEG.  F  ABOVE  ROuM  TEMPERATURE 

BUCKLING  AMPLITUDE  =  .2  15.'  INCHES 

FUNDAMENTAL  FREQUENCY  =  417.5  HZ 


STRES  AT  POINT 
X  Y 


.0  ■■ 

3.0 


6.0 

.  U 


DYNAMIC  STRES 
KSI 

9.371 

8.460 


THERMAL  STRES 
KSI 

—  i6.'83o 

-24.342 


C YALES  TO  FAILURE 

1.32+U5 
1 . 4  0  + 1  b 


50 


-  f"  -  - 

fp- 

* 

TABLE  V 

COMPUTER  PROGRAM  FOR  CALCULATING  ELASTIC  MODULUS 

'JM  ~ 

1 

FUNCTION  ETEMP ( T » IFE) 

2 

C 

3 

c 

THIS  FUNCTION  COMPUTES  ELASTIC  MODULUS  FOR 

??4:i 

4 

c 

ALUMINUM  OR  TITANIUM  ALLOY  AS  A  FUNCTION  OF 

5 

c 

TEMPERATURE 

6 

c 

5- 

7 

c 

T  -  INPUT  TEMPERATURE  -  DEG.  F 

i:-- 

8 

c 

IFF  -  ALLOY  CODE 

9 

c 

=  1  TITANIUM 

10 

c 

=  2  ALUMINUM 

-- 

11 

c 

- 

12 

GO  TO  (100*200) » IFF 

'7- 

13 

r  *************** ******* *********************  ******* 

T5=- 

14 

c 

MATERIAL  6AL'4V  TITANIUM  ANNEALED  SHEET 

IS 

c 

REFERENCE  M IL-HDBK-5B 

16 

c 

TEMPERATURE  LIMITATION  800  DEGREES  F 

.  , 

17 

c 

RT<T<800  F 

if’; 

18 

1Q0 

IF(800-T) 180# 190*150 

-  •, 

19 

190 

ETEMP- (1. 030-0. 000375*T)*16.6E+06 

- - 

20 

RETURN 

I7- 

21 

C 

T>800  F 

-■- 

22 

180 

ETEMP=12.1E+06 

23 

WRITE(6,333) 

24 

RETURN 

?-• 

25 

(3  ************************************************** 

‘^=  -  ~ 

26 

c 

MATERIAL  707S-T6  SHEET 

V, 

27 

c 

REFERENCE  MIL-HDBK-5B 

28 

c 

TEMPERATURE  LIMITATION  600  DEGREES  F 

29 

c 

RT<T<200  F 

V-_- 

30 

200 

1F(200-T)220»210»210 

31 

210 

ETEMP= (1,020-0.00030*T)*10.3E+Q6 

i 

32 

return 

33 

c 

200<T<400  F 

34 

220 

IF(400-T)240#230#230 

tr; 

35 

230 

ETEMP=(0,96-0.OOG70*<T-20O) )*t0.3E+06 

- 

36 

RETURN 

37 

C 

400<T<600  F 

38 

240 

IF(600-T)260#250#250 

39 

250 

ETEMP=(0.82-0.0U16*( r-400) )*10.3E+06 

40 

RETURN 

41 

C 

T>6Q0  F 

42 

260 

ETEMP=0,50*10.3E+06 

43 

WRITE(6,3.33) 

44 

333 

FORMAT ( / #  5x  #  *  UPPER  TEMP  LIMIT  ON  ELAST  MODULUS’ 

45 

1*  EXCEEDED’#/) 

.. 

46 

RETURN 

■*— 

47 

END 

51 


TABLE  VI 


COMPUTER  PROGRAM  FOR  CALCULATING  COEFFICIENT  OF  THERMAL  EXPANSION 


1 

2 

3 

4 

5 

6 
7 
fl 
9 

10 
3  1 
3  2 
13 
lu. 
1  5 
In 

17 

18 

19 

20 
21 
2? 

23 

24 

25 
28 

27 

28 

29 

30 

31 

32 

33 

34 

35 
38 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 
60 

51 

52 

53 

54 


C 

C 

C 

C 

C 

C 

C 

c 

c 

c 


FUNCTION  ALPHA ( T# IFF  > 

THIS  FUNCTION  COMPUTES  COEFFICIENT  OF  THERMAL 

expansion  for  aluminum  or  titanium  alloy  as  a 

FUNCTION  OF  TEMPERATURE 

T  -  INPUT  TEMPERATURE  -  DEG.  F 
IFF  -  ALLOY  CODE 
=1  TITANIUM 
=2  ALUMINUM 


GO  TO  < 100  .200 ) »  IFF 

C  ♦♦♦♦♦♦♦♦♦♦♦t*******'****************  »■**.«**♦***  + 
C  MATERIAL  6AL-4V  TITANIUM  SHErT 

C  ANNEALED 

C  REFERENCE  MIL-HD3K-5B 

C  TEMPERATURE  LIMITATION  1000  DEGRETS  F 

C  RT<T<20C  F 

100  IF(2C0-T>180. 150.150 

150  ALPHA: (4 .45+0. 0C425*T  >  *1 . OE-06 
RETURN 

C  200<T<400  F 

180  IF(400-T) 185# 190. 190 

1  90  ALPHA:  (  4 . 80  +0  .  'inr5*T  )  '*  1  .  OE-06 


RETURN 

C  400<T<1000  F 

185  IF< 1Q00-T) 195. 198.198 
195  WR I TE ( 6 , 500 ) 

500  FORMAT (At 5x. * UPPER  TEMP  LIMIT  ON  COEFF  OF  EXPAn  •» 
1 ’EXCEEDED’ ./) 

198  ALPHA:5.8E-06 


RETURN 

C  *********************************  *************** 

C  MATERIAL  7075-T6  SHEET 

C  REFERENCE  MIL-HDBK-5B 

C  TEMPERATURE  LIMITATION  600  DEGREES  F 

C  RT<T<100  F 

200  IF(10Q-T)220» 210.210 

210  ALPHA: ( 12 . 4+0 . 0050 *T ) *1 . 0 E-06 
RETURN 

C  100<T<300  F 

220  IF(300-T)240. 230.230 

230  ALPHA: (12.9+0.Q0275*(T”100) >*1, OE-06 
RETURN 

C  300<T<400  F 

240  IFI400-T3260.250.250 

250  ALPHA: (13.45+0.00150*( T-300 ) ) *1 .OE-06 
RETURN 


C  T>400  F 

260  ALPHA:13.6E“06 


IF (600-T3280.270.270 
280  WRITEI6.500) 

270  RETURN 
END 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 


TABLE  VII 

COMPUTER  PROGRAM  FOR  CALCULATING  FATIGUE  LIFE 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


SUBROUTINE  SN ( SDYN# STEMP» TEMP* CTF . IFF ) 

THIS  SUBROUTINE  CALCULATES  FATIGUE  LIFE  FOR 
ALUMINUM  OR  TITANIUM  alloy  as  a  FUNCTION  OF 
TEMPERATURE  AND  MEAN  STRESS.  THIS  SUBROUTINE 
IS  BASED  ON  COUPON  FATIGUE  TEST  DATA  AT  ROOM 
AND  ELEVATED  TEMPERATURE. 

ROOM  TEMPERATURE  IS  80  DEG.  F 


SOYN  -  DYNAMIC  STRES  -  KSI  (RMS» 

STEMP  -  THERMAL  (OR  MEAN)  STRESS  -  KSI 
TEMP  -  TEMPERATURE  -  DEG.  F 
CTF  -  LIFE  IN  CYCLES  TO  FAILURE 
IFF  -  ALLOY  CODE 
=  1  TITANIUM 
=  2  aluminum 

GO  TO ( 100  * 200 ) *  IFF 

.***%*****+***************************************** 
MATERIAL  6AL-4V  TITANIUM  SHEET  ANNEALED 
TEMPERATURE  LIMITATION  600  DEG.  F 

100  Cl=12. 58- 0.00376* TEMP 

C2=-5.40+0.00176*TEMP 
ARF=Cl*e2*ALOG10  fSDYN-0. 1*ST£mP) 

CTF=l  0 . **  ARF 
RETURN 

*************************************************** 
MATERIAL  7075-T6  ALUMINUM  SHEET 
TEMPERATURE  LIMITATION  300  DEG.  F 

200  Cl=10.89-0.00584*TEMP 

C2 =-4.89+0. 00 34 7*TEMP 
ARF=C1+C2*ALOG10 (SOYN-O. 1*STEMP) 

ctf=io.*+arf 

RETURN 

END 


I 


C 

c 

c 

c 


53 


TABLE  Vlli 

COMPUTER  PROGRAM  FOR  CALCULATING  SKIN  BUCKLING  TEMPERATURE 


1 

2 

3 

4 

5 
A 

7 

8 
9 

10 
11 
1  2 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 


SUBROUTINE  CTEMP(TCALP»TG. IFF) 

C  THIS  SUBROUTINE  CALCULATES  SKIN  BUCKLING 

C.  TEMPERATURE  FOR  ALUMINUM  OR  TITANIUM  ALLOY 

C  STRUCTURAL  PANELS. 

C 

C  TCALP  -  PRODUCT  OF  BUCKLING  TEMPERATURE 

C  AND  ALPHA  FROM  CALLING  PROGRAM 

C  TC  -  BUCKLING  TEMPERATURE  -  DEG  F  ABOVE 

C  ROOM  TEMPERATURE 

C  IFF  -  ALLOY  CODE 

C  =  4  TITANIUM 

C  =2  ALUMINUM 

C 

TC  =  O.0 

GO  TO  ( 100.200  ) . IFF 

0  ******************%*****%*%**************■********* 
•  C  MATERIAL  6 AL-4 V  TITANIUM  ANNEALED 

C  TEMPERATURE  LlVITA.ION  1000  DEG.  F 

C 

100  C1-4,45E~06 

C2=4 .33E-09 

1=1 

1  C3=60,*C2+C1 

TC=.5*SQRT( <C3/C2)**2*4.*TCALPAC2)-.5*C3/C2 
T=TC+80. 

IF(T-260.)  50.50.2 

2  IF < T-360 • )  3.3.5 

3  Cl=4,9E-06 
C2=2.5E-09 
1*1+1 

IF ( 1-2 )  1.1.50 
5  Cia5.80E-Gfc 
TC=TCALP/C1 
50  TC=TC 
RETURN 
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TABLE  VIII  iCONT) 


36 

37 

38 

39 

40 

41 

42 

43 

44 

45 
45 

47 

48 

49 

50 

51 

52 

53 

54 
56 

56 

57 

58 

59 

60 
61 


Q  ft  | **$,(t**>M:******#  +  *  +  ****SM'*****t1***’*******^:***** 

c  MATERIAL  7075-T6  ALUMINUM  ALLOY 

C  TEMPERATURE  LIMITATION  600  DE6.  F 

c 

200  F1=13.4E*06 
F2X5.0E-09 
11  =  1 

201  F3=80.*F2+F1 

TC=0.50*5QRT<  (P3/F2)**2+>4.0*TCALP/F2>+>0.50*F3/F2 
T=TC+90. 

IFT-100. > 500 » 500*202 

202  lF(T-300. >203,203.204 

203  Fl=l@.625E-06 
F2=0.75E-O9 
11=11+1 

IF ( 1 1-2  >201,201 ,500 

204  JF(T-400. >205.205.206 

£05  Fl=13.0E-06 

F2=1.5E-09 
II* 1 1  +  1 

IF ( 1 1-2  >201 ,201 ,500 
206  Fl=13.6E-06 

TC=TCALP/F1 
500  TC=TC 
RETURN 
END 
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TARLE  IX 

COMPUTER  PROGRAM  FOR  CALCULATING  SECTION  PROPERTIES 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 


SUBROUTINE  PROP ( OPT. B » H. T. A . R J. GAMAT *  X  IP) 

C 

C  SECTION  PROPERTIES 

C  IF  OPT  =  0  ZEE  SECTION 

c  if  opt  a  l  Channel  section 

C  H=  STRINGER  HEIGHT*  CL  TO  CL 

C  B=  FLANGE  WIDTH 

C  T=  STRINGER  THICKNESS 

C  REFERENCE:  AFFDl-TR-71-107 

C  WC  -  WARPING  CONSTANT  ABOUT  SHEAR  CENTER 

C  GAMAT  -  WARDING  CONSTANT  ABOUT  ATTACH  POINT 

C  A  -  CROSS  SECTIONAL  AREA 

C  RJ-  ST.  VENANTS  TORSION  CONSTANT  -  J 

C  XIP  -  POLAR  MOMENT  OF  INERTIA  ABOUT  ATTACH  POINT 

IF ( OPT  1  1*1.2 

C  ZEE  STIFFENER 

1  A=T* (H+2.SB) 

0=H**2*(6.*B+H) 

Ol=T**2*  (3.*Ht-2.*B) 

XXl=(T*t0+01) 1/12. 

02-2 • +B+T 

03=2 .*R-T 

XZI*-<T*H*D2*D3)/8. 

7?I=(T/12. 1*(8.*B**3+H*T**2) 

RJ=<  T**3/3. )*(2.*B+H) 

WC=T*B**3*H**2*{B+2.*H)/'(12.*  <2.*B+H)  ) 

SX=B/2. 

SZ=-H/2. 

D4=SX**2+SZ**2 

D5=D4*A 

X  IPs  »  xx  l  +  2Z  1+05 ) 

GAMAT=WC+ i SZ**2 ) *ZZ 1-2 . *SX*S7*XZ I-H  SX**2 ) *XX I 
RETURN 

C  CHANNEL  SECTION 

2  F=2.*B+H 
XBAR=8**2/F 
F1=6.*B+H 

E=3.*B**2/F1 

CX=E+XBAR 

SXSE+IB/2.I 

SZ=-H/2. 

EX=CX-SX 

A=T*(H+2.*9> 

F2=3 • *H+2 • *B 

XX  I  =  T*  <H**2*F1+T**2*F2  )  /12  . 

F3=12.*H*XBAR**2+8.*B**3 

F4=B-X8AR 

F5=B-2.*XBAP 

Z2  I  =  T* ( F3-24,*XBAR*5wF4  +  12. *B*F5*T+6 . *F4*T*+2  +  T**3 ) /12 . 
RJ=T**3*F/3. 

F6=3,*B+2 ,*H 

WC=T*B**3*H**2*F6/  (  12  .  *F1 ) 

GAMAT=WC+  <SZ**2 ) *Z Z I + ( SX**2 ) *XX I 

F7=(EX**2+SZ**2)*A 

XIP=(XXI+ZZI+F7) 

RETURN 

END 
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VI  -  LIMITATIONS 


Application  of  these  design  procedures  should  be  tempered  with  o  thorough  understanding 
of  their  limitations.  Certain  cf  the  initial  assumptions  stated  during  the  analytical  develop¬ 
ment  were  negated  by  derivation  of  the  empirical  relations.  However,  the  limits  in  the 
physicul  and  environmental  parameters  of  the  experimental  program  then  apply  to  these 
design  criteria,  as  itemized  below. 

A.  Physical  Constraints 

The  bounds  of  the  test  specimen  dimensions  were  used  to  establish  these  limitations. 

These  should  be  treated  only  as  a  guide,  since  the  equations  and  nomographs  are 
normally  valid  beyond  these  limits.  Individual  judgment  mist  be  applied  in  unusual 
cases  where  the  constraints  are  drastically  exceeded,  particularly  in  case  of  the  design 
charts.  The  size  guidelines  are: 

o  Panel  bay  width:  a  =  5  to  9  inches 

o  Panel  bay  aspect  ratio:  b/a  =  1 . 5  to  3 

o  Panel  skin  thickness:  h  =  0.024  to  0.063  inch 

B.  Environmental  Constrai. 


The  acoustic  environment  generally  has  no  restrictions  with  regard  to  the  applicability 
of  the  design  criteria.  Spectrum  levels  below  120  dB  will  normally  result  in  low  dynamic 
stresses  and  a  long  fatigue  life.  The  higher  noise  levels  will  generally  result  in  non¬ 
linear  response,  but  these  effects  are  included  in  the  design  criteria,  since  many  of  the 
test  specimens  exhibited  a  high  degree  of  nonlinearity. 

The  thermal  environment  must  be  nearly  uniform  over  the  surfac;  of  a  panel  boy.  The 
skin  temperature  is  limited  to  the  maximum  temperature  for  which  the  alloys  retain 
significant  structural  properties.  These  limiting  temperatures  are,  for  the  alloys  considered: 

o  7075-T6  aluminum  -  300°F  maximum 

o  6AI-4V  annealed  titanium  -  600°F  maximum 

The  design  life  criteria  are  based  only  on  specific  nominal  temperatures,  requiring  the 
use  of  interpolation  for  intermediate  temperatures,  Extrapolation  beyond  the  temperature 
iimits  ma"  be  permissible  to  some  extent  if  care  is  exercised  and  the  further  temperature 
degradation  effects  are  included. 

Probably  the  most  important  restriction  to  the  design  method  is  in  the  estimation  of  the 
ambient  temperature  and  the  state  of  the  structure  at  this  temperature.  All  thermal  response 
relations  are  referenced  to  the  ambient  temperature  and  the  assumption  that  a  state  of 
stress  equilibrium  exists  (i.e.  ,  no  mean  stresses).  It  is  impractical  at  this  stage  to  give 
guidelines  for  estimating  the  ambient  temperature  state  because  it  is  dependent  on  factors 


such  as  the  length  of  time  ot  a  uniform  temperature  and  external  constraints.  It  should 
be  noted  that  a  change  in  ambient  temperature  over  a  short  interval  constitutes  a  tempera¬ 
ture  change  as  far  as  the  response  relations  are  concerned. 

C.  External  Constraints 

The  external  constraints  imposed  on  the  test  panels  precluded  significant  thermal  expan¬ 
sion  of  the  substructure.  This  is  considered  representative  of  structural  applications  in 
the  direct  flow  path  of  engine  exhausts  or  other  heat  sources,-  where  only  localized  areas 
of  the  structure  are  heated .  The  cri  teria  can  also  be  applied  to  design  applications 
involving  gradual  heating  of  an  entire  structural  area  (i.e.,  supersonic  aircraft),  where  all 
structure  expands  at  about  the  same  rate.  This  corresponds  to  a  relaxation  of  the  external 
constraints  from  those  considered  here.  In  this  case,  the  thermal  buckling  amplitudes 
and  stresses  given  by  the  relations  delineated  herein  will  result  in  a  conservative  design. 

It  should  be  noted  that  the  empirical  results  presented  herein  are  applicable  only  for  the 
case  of  simultaneous  application  of  heat  and  noise.  Alternate  application  of  these  environ¬ 
ments,  wherein  significant  thermal  stress  cycles  are  incurred,  were  not  considered  in  this 
program. 
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APPENDIX  I 


GENERAL  EXPRESSIONS  FOR  STIFFENER  GEOMETRIC  PROPERTIES 
(From  Reference  3) 


The  geometric  parameters  defined  here  are  developed  in  terms  of  o  centroidal  (x,  y,  z) 
coordinate  system.  General  expressions  for  the  cross-sectional  area,  the  area  moments, 
torsion  constant,  and  warping  constants  are  presented  for  zee  and  channel  cross-section 
shapes.  These  parameters  are  defined  as  follows: 

x  -  the  location  of  the  centroid 


e  -  the  location  of  the  shear  center 
A  -  the  cross-sectional  area 


f  2 

/  z  < 

*  *4 

■/> 

■L  =  J" x<t>d  A 


f.A 


-z)2  +  (0,  z  +x)2  , 


(St.  Venont's  Torsion 
Constant) 


>fx  -  L  =  f 


4>26A 


where  <£(x,z)  is  the  warping  function  for  the  cross-section  with  the  pole  taken  as  the 
shear  center. 


NOTE:  The  Symbols  used  in  this  appendix  are  applicable  only  for  this  section.  Since 
they  are  defined  here  or  in  the  following  pages,  they  are  not  inc  luded  in  the  List  of 
Symbols. 


A  =  t(h  +  2b) 


I  =  -j^  [h2(6b  +  h)  +  t2(3h  +  2b)] 


lxz=-^(2b+t)(2b-t) 


I  =~[8b3+ht2] 
zz  I  z 


J  =  y[2b+h) 

^  _  tb3h2{b  +  2h) 
1  e  12(2b  +  h)_ 


FIGURE  I  -  1  GEOMETRIC  PROPERTIES  -  ZEE  SECTION 
(FROM  REFERENCE  3) 
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x  =  b2/(2b  +  h) 

e  =  3b2/(6b  +  h) 

A  =  t(h  +  2b) 

lKX  =  ^CH2(6b+h)  +  t2(3h+2b)] 
I  =  0 


-  I  =  J2  [12h^2  +  8b3  -  24xb(b  -  x)  +  12b(b  -  2*)t  +  6(b  -  ^)t2  +  t3] 

J  =  ~  t3(2b  +  h) 

tb3h2(3b  +  2h) 

1 «  12(6b  +  h) 


NOTE:  For  worping  constants,  the  pole  is  token  at  the  shear  center. 
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FIGURE  I  -  2  GEOMETRIC  PROPERTIES  -  CHANNEL  SECTION 
(FROM  REFERENCE  3) 
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